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A crucial part of aircraft design is ensuring that the required margin for flutter is 
satisfied. A trustworthy flutter analysis, which begins by possessing an accurate dynamics 
model, is necessary for this task. Traditionally, a model was updated manually by fine tuning 
specific stiffness parameters until the analytical results matched test data. This is a time 
consuming iterative process. NASA Dryden Flight Research Center has developed a mode 
matching code to execute this process in a more efficient manner. Recently, this code was 
implemented in the F-15B/Quiet Spike™ (Gulfstream Aerospace Corporation, Savannah, 
Georgia) model update. A build-up approach requiring several ground vibration test 
configurations and a series of model updates was implemented in order to determine the 
connection stiffness between aircraft and test article. The mode matching code successfully 
updated various models for the F-15B/Quiet Spike™ project to within 1 percent error in 
frequency and the modal assurance criteria values ranged from 88.51-99.42 percent. 
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moment of inertia about the y-axis with respect to the z-axis, lb-in 2 
moment of inertia about the y-axis with respect to the z-axis, lb-in 2 , target 
moment of inertia about the z-axis with respect to the x-axis, lb-in 2 
moment of inertia about the z-axis with respect to the x-axis, lb-in 2 , target 
second index term 
objective function 
stiffness matrix 
mass matrix 

Modal Assurance Criteria 
number of degrees of freedom 
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I. Introduction 

W hen a new aircraft is designed or a modification is made to an existing aircraft, the aeroelastic properties of 
the aircraft need to be examined to ensure it is flightworthy. Evaluating the aeroelastic properties of a new or 
modified aircraft can include performing a variety of analyses, such as modal and flutter analyses. To produce 
accurate results from these analyses, it is imperative to work with finite element (FE) models that have been 
validated by or correlated to ground vibration test (GVT) data. Updating an analytical model using measured data is 
a challenge in the area of structural dynamics. The analytical model update process encompasses a series of 
optimizations that match analytical frequencies and mode shapes to the measured modal characteristics of a 
structure. In the past, the method used to update a model to match testing data has been ‘trial and error.’ This is an 
inefficient method; running a modal analysis, comparing the analytical results to the GVT data, manually modifying 
one or more structural parameters (mass, center of gravity (CG), inertia, area, etc.), rerunning the analysis, and 
comparing the new analytical modal characteristics to the GVT modal data make up a time-consuming process. If 
the match is close enough (close enough defined by the analyst’s updating requirements), then the updating process 
is finished. If the match does not meet updating requirements, then the parameters are changed again and the process 
is repeated. Clearly, this manual optimization process is highly inefficient for large FE models and/or a large number 
of structural parameters. 

The NASA Dryden Flight Research Center (Edwards, California) (DFRC) has developed, in-house, a mode 
matching code that automates the above-mentioned optimization process. Moreover, it facilitates the development of 
reduced order models, which simplify analyses and therefore, reduce computation time. The DFRC in-house mode 
matching code reads mode shapes and frequencies that were measured from a GVT to create the target model. It also 
reads the current analytical model as the ‘initial guess,’ as well as the variable parameters and their upper and lower 
limits as assigned by the analyst. It performs a modal analysis on this initial model and modifies its variable 
properties to create an updated model that has similar mode shapes and frequencies to those of the target model. The 
mode matching code outputs frequencies and modal assurance criteria (MAC) values that allow for the quantified 
comparison of the updated model versus the target model. 

The mode matching code assists in determining the airworthiness of the aircraft flown at DFRC. Recently, this 
code was implemented on a project flown on the DFRC F-15B airplane supersonic flight-testing platform, the 
Gulfstream Aerospace Corporation (Savannah, Georgia) Quiet Spike™ (QS). This code allowed for a build-up 
approach that was necessary because the flight test article itself was unavailable and allowed for progress to be made 
in the analyses even without the flight hardware. 

The objective of this work was to use the mode matching code to determine the connection stiffness on the 
aircraft radar bulkhead where the QS and the F-15B airplane mate, such that when the FE model of the flight QS 
was validated, it could be readily mated to the aircraft FE model. This approach was established to enable the project 
to progress the model development and flutter predictions without the actual QS flight test hardware. 


II. Mathematical Background 

Three optimization phases are used to refine a model with the DFRC mode matching code. The mass properties 
are modified to match measured mass properties, the mass matrix is orthogonalized, and then the natural frequencies 
and mode shapes are matched. Design variables for the optimization can include structural sizing information 
(thickness, cross-sectional area, area moment of inertia, torsional constant, etc.), point properties (lumped mass, 
spring constant, etc.), and material properties (Young’s modulus, etc.). 
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Phase 1 : Mass properties 

Matching mass properties by using lumped masses (CONM2) 1 reduces the computational effort 
required to match stiffness properties and ensures that a physically meaningful and accurate solution is 
achieved. In this first optimization phase, 10 objective functions are minimized to match analytical to 
measured mass properties. The mass properties that are updated include total mass, center of gravity 
location, and mass moments of inertia. As an objective function is minimized, a constraint function is 
implemented to prevent the updated variable from changing in the minimization of the following objective 
function. The constraint functions limit the amount of change allowed in the previously optimized 
variables. The objective functions and constraints applied in this phase are listed in order in Table 1. This 
optimization phase places the model within the vicinity of a physically viable solution (feasible region), 
facilitating the second and third optimization phases. 

Table 1. Optimization Problem Statement for Phase 1 of Model Update. 


Statement 

Objective Function 

Constraint Function 

1 

J ! = (W-W G ) 2 /W G 2 

Unconstrained 

2 

J 2 = (X-X g ) 2 /X g 2 

w 

< e 

3 

J 3 = (Y-Y g ) 2 /Y g 2 

W < e 

i=l,2 

4 

Ja = (Z-Z g ) 2 /Z g 2 

|Ji|<e i 

= 1, .. 3 

5 

Ji = (Ixx _ Ixxg)”/Ixxg~ 

|/i|<e i 

= 1, ... 4 

6 

J (, = (Iyy-Iyyg) /Iyyg 

|/i|<e i 

= 1, ... 5 

7 

Ji~ (Izz-Izzg)”/Izzg _ 

|Ji|<e i 

= 1, ... 6 

8 

J» = (Ixy-Ixyg) /Ixyg 

Wl<£ i 

= 1, ... 7 

9 

J 9 = (Iyz-Iyzg)~/Iyzg 

1-Til < £ i 

= 1, ... 8 

10 

Jio ~ (Izx-Izxg)~/Izxg 

|«Ti| < £ i 

= 1, ... 9 


Phase 2: Mass Matrix Orthogonalization 

The orthonomialized mass matrix is shown in Eq. (1). 

M = O g T T T MT<I> g (1) 

The off-diagonal terms of this orthonomialized matrix are reduced to improve mass orthogonality. In 
Eq. (1), the mass matrix, denoted as M, represents the analytical mass matrix, while 3> G represents the 
target Eigen matrix that is obtained from the GVT. The target Eigen matrix, <£> G , remains unchanged during 
the optimization. The transformation matrix, T, is based on Guyan reduction 2 or improved reduced system. 3 
The optimization problem statement then becomes: 

Minimize 


Jjj = jr (Mi /) 2 (Off-diagonal terms of M) (2) 

Such that |J ; | < s, where i = 1, 2... 10(7; values are given in Table 1), and n denotes the number of modes 
that will be correlated in Eq. (2). The positive-definiteness characteristic of the lumped mass elements also 
acts as a constraint during this optimization procedure. 

Phase 3: Frequencies and Mode Shapes 

In this phase, certain modes can be given a higher level of importance by applying a scaling factor, (5, in 
Eqs. (3) and (5). This scaling factor helps to more closely tune the frequencies of these analytical modes of 
interest after the mode shape vectors have been matched to within required accuracy. 
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Two optimization methods are used in this phase. In the first method, the objective function considered 
combines an index which identifies normalized errors between test (CO;) and analytical (Q;) frequencies with 
a second index that defines the total error associated with the off-diagonal terms of the orthonormalized 
stiffness matrix. The optimization problem statement then becomes: 

Minimize 


n 

J\2 = X 


i=l 


/„ \ 2 

Qi~ coj 

V / 


i (& ,) 2 


( 3 ) 


Such that |Jj < e for i = 1, 2... 11. The matrix K is a result of the matrix multiplication shown in Eq. (4), 
where the stiffness matrix, K, is the analytical stiffness matrix. 

K =0 (j r T KT0 (| (4) 

In the second optimization method, the objective function combines the same frequency error index used 
above, with a second index which defines the total error between test and analytical mode shapes at 
specified sensor locations. The optimization problem statement then becomes: 

Minimize 


J [ 2 ~ X 


' Qi~cPi x 


i=l v 


+ aX(0>i-OiG)“ 

i=l 


( 5 ) 


Such that \Ji\ < e for i = 1, 2... 11, and m denotes the number of degrees of freedom measured. The 
frequencies CO; and mode shapes <t>j are calculated based on the direct-iterative Eigen solution technique for 
large eigenvalue problems 4 : 


1) Compute K t and M T : 

2) Compute F: 

3) Compute R: 

4) Compute K R and M R : 

5) Compute ( 0 ; and 'E : 

6) Compute <& : 


K t = T t KT and M x = T t MT 
F — MjOq 

f = k t r 

Kr = RK X R and M R = RM X R 
Kr'E = co i 2 M R 'F 

<& = R»P 


For the direct-iterative Eigen-solution technique, the analytical stiffness and mass matrices, K and M, 
are correlated, while the measured mode shapes, <J> G , remain unchanged during the model updating process. 
Notice that the transformation matrix, T, in Eqs. (1) and (4) is calculated using the analytical mass and 
stiffness matrices M and K. Therefore, the constraint equation |/u| < e is needed for the third optimization 
phase. 

Current model updating approaches consist of updating only stiffness properties and assume that the 
modeled mass distribution is accurate. This can create erroneous stiffness properties in a model. The DFRC 
mode matching code is capable of updating both mass and stiffness in a model by updating mass properties 
such as total mass, CG, and mass moments of inertia, improving the orthogonalization of the mass matrix, 
and matching frequencies and mode shapes to test data. Updating both mass and stiffness properties of a 
model will provide a higher degree of confidence in a model’s resemblance of the actual hardware. 


III. Application 

An example of an application of the mode matching code is the F-15B research test bed. The DFRC possesses a 
modified F-15B airplane that is used as a test bed for supersonic flight experiments. Traditionally, the FE model of 
the flight test article is generated and updated by GVT data and then mated to the updated F-15B FE model. A GVT 
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is then conducted with the test article mated to the aircraft. The mated FE model is correlated to this mated GVT 
data with the going-in assumption that, at this point, the only unknown variable is the interface stiffness between the 
test article and the aircraft. 

Recently, Gulfstream Aerospace Corporation (GAC) desired to test a modified nose, the Quiet Spike™, on the 
DFRC F-15B 5 airplane. Because of flight test article availability and a compressed project schedule, a non- 
traditional testing and model update method had to be implemented. A build-up approach was devised such that a 
Mock QS would be designed and fabricated to replace the flight QS in the tests and analyses. A detailed FE model 
of this Mock QS would be created. A reduced order model of this detailed FE model would be created in the form of 
a beam model to expedite the model update process. The beam model would then be matched to data provided by 
the first GVT using a strongback. This beam model would be mated to the F-15B model and the interface stiffness 
would be updated by a second GVT using the airplane. The flight QS would be tested and its model would be 
updated separately at GAC. It would then be mated to the aircraft model and its modal analysis results would be 
compared to frequencies and mode shapes from a final GVT of the flight QS mated to the F-15B airplane. 
Ultimately, the goal of this build-up approach was to completely eliminate the need for a model update of this final 
flight test configuration. 

The flight QS measures 30 ft in length when fully extended and 20 ft when fully retracted. Figure 1 shows the 
GAC flight QS. 



Radome 

bulkhead 



Figure 1. Flight QS Description. 

The Mock QS was designed with the purpose of replacing the flight QS in the build-up GVT approach. For the 
Mock QS to fulfill this purpose, it was necessary that its design meet the following requirements: 1) similar weight 
to flight QS weight, 2) similar CG to flight QS CG, and 3) similar Ixx for the third segment of the flight QS. A final 
design requirement imposed on the Mock QS was that it was to interface with the strongback and F-15B radar 
bulkhead in the same manner and with the same hardware as the flight QS installation assembly. A study was done, 
which considered different materials, beam types, thicknesses, diameters, and lengths to determine the best design 
for the Mock QS. The Mock QS measures 19 ft and is shown in Fig. 2. 
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Figure 2. Mock QS Description. 

Several model updates coupled with a series of four GVTs were performed to acquire a correlated mated 
F-15B/QS FE model that could be used in flutter analyses. Figure 3 shows the various models of the Mock QS that 
were developed to accomplish this task. The first FE model created was a detailed FE model of the Mock QS, which 
was created in Patran and MSC.Nastran (Santa Ana, California) using the Mock QS computer aided design (CAD) 
models. It consisted of beam, plate, and solid elements and was fixed in a manner that represented how the flight QS 
mounted to the airplane. A modal analysis was performed on the detailed FE model. The results from this modal 
analysis rendered the target frequencies and mode shapes for the beam model, which is shown on the far right in 
Fig. 3. 



Figure 3. Mock QS model development. 

The original beam model was developed in Patran and MSC.Nastran. The lumped mass and beam elements were 
created manually based on the numerical mass properties of the detailed model. Because the model update process 
requires much iteration, using a reduced order model decreases the computational effort required to match a model’s 
structural properties to test data. This beam model decreased the computational effort in all analyses including the 
flutter analyses required at the conclusion of the model update process. Furthermore, matching the beam model to 
the detailed model first provided a good ‘initial guess’ for the Mock QS beam model and would, therefore, facilitate 
correlating the beam model to test data. 

The next step in this task was to correlate the Mock QS beam model to the actual Mock QS structure. This was 
accomplished in the first GVT, which consisted of attaching the Mock QS boom to a strongback and exciting it in 
the lateral and vertical directions as shown in Fig. 4. The results from this GVT were used to update only the 
structural properties of the beam model. The interface stiffness was not updated until the second GVT. Frequencies 
and mode shapes were extracted from all the following GVT’s to provide target data for model validation. 6 
Attaching the Mock QS structure to the Strongback provided a nearly rigid foundation for the boom that ensured the 
boom modes were excited well. This also matched well with the boundary conditioned modeled. 
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Figure 4. Mock QS ground vibration testing on the Strongback structure. 


When mating the updated Mock QS beam model to the updated aircraft beam model, it was assumed that any 
discrepancies between the numerical and the test mode shapes and frequencies could be attributed to the unknown 
connection stiffness between the test article and aircraft. The second GVT would empirically provide this 
connection stiffness. The aircraft was constrained laterally and vertically, such that the aircraft modes would be 
decoupled from the Mock QS modes. The Mock QS was again excited both laterally and vertically. This GVT 
configuration is illustrated in Fig. 5. 



Figure 5. Mock QS boom ground vibration testing on the F-15B airplane. 


The third GVT consisted of the flight QS mounted on a strongback that provided the modal parameters to update 
the flight QS model. The fourth and final GVT served as a final check for the mated F-15B/QS frequencies and 
mode shapes. 
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IV. Results and Discussion 

Because of the simplicity of the structure and the fidelity of the detailed FE model of the Mock QS, there was 
confidence in its representation of the actual Mock QS hardware. Therefore, in performing the first model update, 
the target modal values were provided by the numerical results of the detailed model to update the beam model. The 
first three modes of the detailed model or, target modes, are shown in Fig. 6. 




Mode 1 - Lateral bending, 5.675 Hz 


Mode 2 - Vertical bending, 6.602 Hz 



Mode 3 - Torsion, 7.876 Hz 


Figure 6. Detailed Mock QS inode shapes and frequencies. 


The first three modes of the detailed Mock QS were lateral bending at 5.68 Hz, vertical bending at 6.6 Hz, and 
torsion at 7.88 Hz. Table 2 quantifies the similarity in the modal parameters between the beam and detailed Mock 
QS models before the beam model was updated. 

Table 2. Comparison of Frequencies and MAC Values before Update. 


Mode 

Detailed FE Model (Hz) 

Beam Model (Hz) 

MAC (%) 

Freq. Error (%) 

1 

5.675 

7.285 (3 rd detailed mode) 

0.05 

28.37 

2 

6.602 

7.515 (1 st detailed mode) 

0.02 

13.83 

3 

7.876 

8.977 (2 nd detailed mode) 

0.03 

13.98 


Table 2 indicates that, initially, the beam model did not correlate well with the detailed model. The beam 
model’s modes were ordered differently - 1) torsion at 7.285 Hz, 2) lateral bending at 7.515 Hz, and 3) vertical 
bending at 8.977 Hz. Consequently, this resulted in very low MAC values and high errors in frequency. The higher- 
than-expected frequencies imply that either the masses were underestimated or the stiffnesses were overestimated in 
the beam model. 

As discussed in the previous section, the first phase of updating the beam model consisted of updating the mass 
properties. The optimization variables were assigned upper and lower limits of + 10 percent of the original values 
and included the mass values, CGs and inertias of the lumped mass elements. All mass properties were matched to 
within 0.5 percent error with the exception of the y-CG because the detailed model possesses a y-CG at zero. A 
comparison for all mass properties is shown in Table 3. 

Table 3. Comparison of Mass Properties after Update. 



Mass, 

lb 

x-CG, 

in 

y-CG, 

in 

z-CG, 

in 

fxx, 

lb-in 2 

Iyy, 

lb-in 2 

Izzj 

lb-in 2 

1 YZ, 

lb-in 2 

l/x, 

lb-in 2 

Detailed 

400.01 

114.67 

0.00 

108.08 

2.60E+04 

2.16E+06 

2.15E+06 

7.40 

2.41E+04 

Beam 

400.59 

114.55 

0.03 

108.10 

2.60E+04 

2.16E+06 

2.15E+06 

7.38 

2.41E+04 

Error 

-0.14% 

0.10% 

NaN 

-0.02% 

-0.14% 

0.01% 

0.01% 

0.21% 

0.10% 
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The second phase of this model update consisted of producing an orthogonalized mass matrix. The variables 
were the same as in the previous phase, however, with the implementation of the constraint functions, the variables 
were limited to changing in such a way that the objective functions errors from the previous phase would not 
increase to a value greater than e (See Table 1). At the conclusion of this phase, the off-diagonal term in the mass 
matrix of the greatest value was 0.0135. The resulting orthogonalized mass matrix was: 

1.0000 -0.201 IE-03 0.1352E-01 

-0.201 IE-03 1.0000 0.3827E-02 

0.1352E-01 0.3827E-02 1.0000 

The third phase of implementing DFRC’s mode matching code matches the beam model frequencies and mode 
shapes to the detailed model results. For this phase, the design variables were expanded to include element stiffness 
properties such as beam cross sectional area, area moment of inertias, and torsional constants. This set of design 
variables was also allowed to range + 10 percent to complete the model update procedure. Table 4 demonstrates 
how the DFRC in-house mode matching code correlated the beam model to the detailed model. 

Table 4. Comparisons of Frequencies after the First Update. 


Mode 

Detailed FE Model (Hz) 

Beam Model (Hz) 

MAC (%) 

Freq. Error (%) 

1 

5.675 

5.674 

98.29 

-0.02 

2 

6.602 

6.607 

98.14 

0.08 

3 

7.876 

7.867 

78.64 

-0.11 


The updated beam model showed good correlation to the detailed model. The first two mode shape vectors were 
matched to less than two percent error, and the frequencies for all three modes were matched to within one percent 
error. The third mode shape, which was the torsion mode, did not possess as good a MAC value as the first two. 
Because of the geometry of the QS and the additional torsional stiffness provided by the radome, the first torsion 
mode of the actual QS does not appear until about 100 Hz. Therefore, it was deemed not critical to match the torsion 
mode with as much rigor as the lateral and bending modes, and a MAC value of 78 percent sufficed for this mode. 

Figure 7 shows the first three mode shapes of the Mock QS beam model after being updated. 



Figure 7. Mock QS beam model mode shapes and frequencies. 


The beam model was next updated to match test data that resulted from the first GVT. This first Mock QS GVT 
required that the Mock QS be mounted to a strongback. The first three modes of the Mock QS were 1) lateral 
bending at 5.583 Hz, 2) vertical bending at 7.752 Hz, and 3) torsion at 8.105 Hz. Updating the beam model to the 
detailed model first, made updating the beam model to test data faster and less computationally intensive. Table 5 
compares the beam model to the strongback GVT results after only updating the mass properties. In this case, the 
only mass properties that were measured were total weight, and x-CG. The numerical values were kept for all other 
mass properties. 


9 

American Institute of Aeronautics and Astronautics 


Table 5. Comparison of Frequency and MAC Values at Update Start. 


Mode 

GVT (Hz) 

Beam Model (Hz) 

MAC (%) 

Freq. Error (%) 

1 

5.583 

5.643 

99.29 

1.08 

2 

7.752 

7.529 

99.49 

-2.88 

3 

8.105 

8.142 

67.66 

0.45 


Because the beam model had already been through one update, its modal parameters already correlated well with 
the measured modal parameters. The design variables were treated similarly as before. They were assigned upper 
and lower limits of ± 10 percent of the original values and included the mass values, CG’s and inertias of the 
lumped mass elements. The total weight and x-CG were updated to within a one percent error. 

The second phase of the model update, mass matrix orthogonalization, was also treated as previously done. The 
design variables were the same as in the previous phase. At the conclusion of this phase, the off-diagonal term in the 
orthogonalized mass matrix of the greatest value was 0.0737. The resulting orthogonalized matrix was: 

1.0000 -0.1537E-01 0.7368E-01 

-0.1537E-01 1.0000 -0.6174E-01 

0.7368E-01 -0.6174E-01 1.0000 

Finally, the beam model was updated to have frequencies and mode shapes that matched the measured data. The 
design variables were again expanded to include element stiffness properties such as beam cross sectional areas, area 
moment of inertias, and torsional constants that were allowed to range ± 10 percent. Table 6 shows how the DFRC 
in-house mode matching code correlated the beam model to the actual Mock QS boom. 

Table 6. Comparisons of Frequencies after Update to GVT Data. 


Mode 

GVT (Hz) 

FE Model (Hz) 

MAC (%) 

Error (%) 

1 

5.583 

5.594 

99.25 

0.20 

2 

7.752 

7.764 

99.42 

0.15 

3 

8.105 

8.093 

88.51 

-0.15 


The resultant numerical modal parameters showed an excellent correlation to the measured modal parameters for 
the Mock QS test article. The first two mode shape vectors were matched to less than one percent error, and the 
frequencies for all three modes were matched to within 0.5 percent error. Regarding the third mode shape, we see 
the same results as before. This mode, the torsion mode, did not possess as good a MAC value as the first two. 
Again, this mode was regarded as not being of great importance. 

The final model update required that only the spring element stiffnesses that connected the Mock QS hardware to 
the aircraft be updated. Initially, they were assumed to be nearly rigid and therefore, were assigned high values. The 
second GVT would provide their true values. The second GVT consisted of mating the Mock QS boom to the F-15B 
airplane. The airplane was constrained vertically by supporting it with aircraft jacks and laterally by use of an 
erector set. This configuration decoupled the aircraft modes from the boom modes. The measured modes and 
frequencies of this configuration were lateral at 6.06 FIz and vertical at 7.75 Hz. The torsion mode was not extracted 
for this configuration. Table 7 computes the frequency errors between the numerical and test modes. 


Table 7. Comparison of Test and Numerical Modes of the Mated Mock QS/F-15B Airplane. 


Mode 

GVT (Hz) 

FE Model (Hz) 

Freq. Error (%) 

1 

6.064 

5.583 

8.62 

2 

7.756 

7.752 

0.05 
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The first airplane test mode was at a higher frequency than the first numerical mode by almost 0.5 Hz. This was 
not expected, because the mated model already assumed a nearly rigid connection between the Mock QS and the 
airplane. The higher frequency of the first airplane test mode indicates that the airplane was stiffer laterally than the 
strongback, which can be explained by the high inertia of aircraft compared to the strongback and the erector set 
providing a true lateral constraint. This also indicates that the assumption of a nearly rigid connection was accurate 
and is further supported by the fact that the second numerical mode frequency already matched the second test mode 
frequency to within a tenth of a percent. 

After the flight QS model was updated by test data, it was mated to the F-15B model using this nearly rigid 
connection stiffness. The project then proceeded to perform flutter analyses using the mated model. This final mated 
model matched well to GVT data, and therefore, no model update was necessary, and the flutter analysis did not 
need to be redone. The QS team could quickly proceed to the flight testing portion of this project. 

V. Conclusions 

The NASA DFRC has developed an in-house mode matching code that updates both the mass and stiffness 
properties of a model. This method renders a more accurate model than the method of only updating stiffness 
properties. The three phases of optimization that the code employs are 1) updating mass properties, 
2) orthogonalization of mass matrix, and 3) updating frequency and mode shape. The code has been applied in the 
Gulfstream Aerospace Corporation QS project flown on the DFRC F-15B supersonic flight test bed. The build-up 
approach used to update the QS/F-15B mated model was completed using this code. During this project, it was used 
to perform two model updates. It successfully updated each model used for each mode of interest to within one 
percent error. This code was applied in such a way that allowed the project to proceed with flutter analyses, such 
that when the flight test article arrived at Dryden, the project could quickly move into flight testing. 
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